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Abstract. We investigate magnetic properties of strongly interacting 
bosonic mixtures confined in one dimensional geometries, focusing on 
recently realized 87 Rb- 41 K gases with tunable interspecies interactions. 
By combining analytical perturbation theory results with density-ma- 
trix-renormalization group calculations, we provide quantitative esti- 
mates of the ground state phase diagram as a function of the rele- 
vant microscopic quantities, identifying the more favorable experimen- 
tal regimes in order to access the various magnetic phases. Finally, we 
qualitatively discuss the observability of such phases in realistic setups 
when finite temperature effects have to be considered. 



1 Introduction 



Experimental advances in the preparation and manipulation of cold atoms and mole- 
cules trapped in optical lattices [l][2] have revived the theoretical interest towards 
models and problems that in the past played a fundamental role for the description 
of possible new states of matter, but lacked physical realization in solid state setups. 
Indeed, thanks to the possibility of using both fermionic and bosonic type of atoms 
and to the high tunability of the interaction shape and parameters, the phenomena 
that one can access with these systems are the most varied, ranging from standard 
superfluidity to BEC-BCS cross-over, to Mott physics and dynamical processes [2]-[5]. 

Among the various advantages in dealing with cold atoms, it is worth men- 
tioning the possibility of trapping them in highly anisotropic optical lattices, thus 
realizing systems with different geometries and dimensions. In particular, strongly 
anisotropic lattices allow to realize systems in which the atoms are forced to live in 
one-dimensional (ID) lattices, i.e. on a chain |2)m. 

From a theoretical point of view, physics in ID is of particular interest not only 
because it may realize toy models for higher dimensional problems, but also because 
some physical effects, such as quantum ones, are much stronger and give rise to 
new phenomena and new states of matter [Hp]- Analitically, one can use powerful 
low-energy field theories to describe fermionic, Dosonic, and spin models in terms of 
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collective bosonic degrees of freedom [7||8j. Also, under suitable hypothesis, one can 
establish an equivalence between fcrmionic and spin degrees of freedom. 

Theoreticians have long been studying models by switching from a bosonic to a 
fermionic representation, from a fermion to a spin one, or viceversa. These mappings 
allow to see how new phases of matter may stabilize for certain values of the coeffi- 
cients appearing in the Hamiltonian. This is the case, for example, of the emergence 
of magnetic ordered phases in fermionic and in bosonic models, a problem that was 
studied long ago in the seminal work by Mott. It was first demonstrated that an in- 
sulating antiferromagnetic phase may arise in an apparently simple fermionic model 
such as the Hubbard one (see [9 10 for reviews). The analysis has then been ex- 
tended to more complicated systems, such as mixtures of fermionic species described 
by the two-species Hubbard model, where also more exotic phenomena may appear 
(singlet superconductivity, FFLO phase, super-counter- flow, etc.) (see Refs. [2||6 11 
for a complete review). 

As for the bosonic case, it was first in Refs 



12 13 that it was shown that the 



Bose-Hubbard model could admit a quantum phase transition from superfluid to an 
insulating magnetic Mott-like phase. Since then, much attention has been devoted to 
understand the features of this transition in arbitrary dimensions or in more complex 
models describing for example mixtures of bosonic species 



14-24 



In recent years, these theoretical predictions have become accessible to experimen- 
tal checks in cold atoms experiments. Starting from the breakthrough demonstration 
of Mott-insulator to superfluid transition in a gas of Rb atoms 25 , strongly correlated 



regimes have been systematically accessed in ultr acold g as experiments, promi nen t 
examples being the realization of Tonks-Girardeau 26p7 and super- Tonks gases 



the observation of ID Mott 29 



and Berezinskii-Kosterlitz-Thouless transitions 
[211311. 



and various dynamical processes 

Even more intriguing is the possibility of investigating magnetism in such control- 
lable setups, as recently shown in the case of the Ising model in Ref. [32]. In partic- 
ular, two-species bosonic and fermionic gases represent ideal setups for the search of 
strongly correlated magnetic states of matter, since once the charge degrees of free- 
dom are gapped in a Mott insulating region, the many-body dynamics is dominated 
by pseudo-spin degrees of freedom 14 15 17 , as in the case of the aforementioned 



Hubbard model 10 



In this work, motivated by impressive experimental achievements in tuning and 
controlling low-dimensional heteronuclear bosonic mixtures |33[|34| , we present a sys- 
tematic and quantitative investigation for the realization of magnetic phases in a 
feasible setup of a Bose-Bose mixture of cold atoms trapped in one-dimensional op- 
tical lattices. In Sec. [2j we will describe how to model the system via a two-species 
Bose-Hubbard Hamiltonian, whose effective parameters are determined by micro- 
scopic properties, such as the depth of the lattice potential and the Feshbach reso- 
nance scattering length. We will discuss this relationship in details for Rb-K mixtures 
which, thanks to their high tunability 33 , allow to span a wide range for both the 



hopping and the interaction coefficients. 

In Sect. 3, we will discuss the theoretical background to study the emergence of 
magnetic phases in such a system. To describe the physics deep inside the insulating 
Mott regions, where the single site density is fixed to integer values, we will study 
the two species Bose-Hubbard Hamiltonian in the strong coupling regime via a per- 
turbative scheme which is a generalization of the Schrieffer- Wolff transformation [35] . 
At integer fillings, the so obtained effective Hamiltonian can then be mapped onto 
a spin model, which turns out to be the spin-1/2 XXZ-cham for filling one particle 
per site and the A — D spin-1 Hamiltonian for filling two. Both these models are 
paradigmatic for the description of many body systems that display quantum phase 
transitions [8|k37 38 . As it is well known, for strong anisotropies, the former one 
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admits a ferromagnetic and an antiferromagnetic regime. The latter has a very in- 
teresting zero-temperature phase diagram, with different critical and massive phases, 
which include ferromagnetic, Ising-like as well as spin-0 ground states. In particu- 
lar, it displays a Haldane phase, which is determined by the breaking of a non-local 
symmetry and described by non-local order parameters, the so-called string-order 
ones. In both cases, we will address how feasible is to reach these regimes from an 
experimental point of view. 

In Sec. [4] we will study the Bose-Hubbard Hamiltonian in the spin-1/2 magnetic 
sector by means of numerical simulations based on the density-matrix-renormalization- 
group (DMRG) algorithm [39|40| . Our calculations allow to determine a quantitative 
phase diagram for the experimentally relevant case of Rb-K mixtures. On the basis of 
this analysis, we will be able to predict the emergence of the various massive phases as 
the parameters change in an experimentally feasible range, and estimate the critical 
temperatures needed to achieve such phases in experiments. 

Finally, in Sec. 5, we draw our conclusions and discuss possible developments. 



2 Model Hamiltonian and experimental regimes 

Two-species Bose mixtures strongly confined by optical lattice potentials in two- 
directions and subject to an additional periodic potential along the third direction 
are well described by the following microscopic Hamiltonian: 



h 2 d 2 

-—^+V(x)+p a (x) 



ip a (x)+ 



01=1,2 

+ / dxdx' p a {x)p a (x')U a 5{x — x')+ (1) 

J a 

+ J dxdx' pi{x)p2{x')U\25{x — x') 

Here, ^ a ,Tp a are bosonic creation/anihilation operators, m a are particle masses, 
V(x) = Fsin 2 (27rcc/A) is the optical lattice potential along the wire of depth V and 
wavelength A, and p{x) the confining potential along the x-direction. The intra-species 
interaction strengths U a = 2huj± ta aiD. a are related to the intraspecies ID scattering 
length a\D,ai with ujj_. a = omtP \Z^±/^r,a being the transverse confinement fre- 
quency depending on both the optical lattice spacing d = A/2 and the intensity of the 
transverse field V± in units of the recoil energy E r ^ a = Sr £ rf2 . Finally, the inter-species 
interaction U12 is usually tunable via Feshbach resonances [41] , and is proportional 
to the inter-species scattering length axon- In the limit of sufficiently deep lattice 
potentials along the wire direction, V/E r , a > 6, the system is well modeled by an 
effective Bose Hubbard Hamiltonian: 

Hbb = - ^2 t a y^(at g a, + i |g + h.c.) + U ab ^ n ita n^ b 

<t i i 

+ n iA n i,<? - l ) + Y ^ Y ~ L / 2 ) 2 

(Ti a i 

where the first term represents tunneling processes of both species a = 1,2, the second 
and third are the inter- and intra-species interaction, and the last one is a species- 
dependent confining potential. The operators at ,a,j t<x satisfy bosonic commutation 
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relations, and n^ a = a\ G a,i^ a - The effective Hubbard parameters can be derived from 
microscopic quantities by expanding the single particle wave-functions in the Wannier 
basis, thus obtaining the following relations [21: 



t„ = 



AE r 



V 

E r .n 



3/4 



exp 



-2 



V 

E r .n 



1/2' 



u„ = 



A 



V 

E rn 



1/4 



(2) 
(3) 



Once a specific experimental setup is considered, the relevant Hamiltonian parame- 
ters can be directly shaped by means of external fields as follows: the hopping rate 
ratio is tuned by considering, e.g., species dependent lattices, or, in case of species 
independent one, by just changing the lattice depth (which affects more drastically 
heavier particles). Interactions can be then tuned by Feshbach resonances; in par- 
ticular, leaving intra-species interactions unmodified, it is possible to span a broad 
range of interaction strengths by means of an external magnetic field, or, eventually, 
employing confinement induced resonances. 



2.1 Experimental parameters for Rb-K mixtures 

As a case study, we will focus from here onwards on 87 Rb- 41 K mixtures. This choice 
is justified, from the one hand, by the notable experimental achievements already 



demonstrated with these mixtures such as, e.g., full tunability of interactions 33 
and, on the other hand, by the relative large mass-imbalance, which allows to span 
parameters regimes ranging from intermediate to large hopping ratios. Here, we con- 
sider optical lattices with wavelength A/2 = 532nm (as reported in [33]), and define 
dimcnsionless optical lattice depths s = V/E Tt K and sj_ = V_i_/E r ,K along the x and 
y — z directions respectively. 

Typical experimental parameters obtained from Eq. [2] [3] are presented in Fig. [I] 
The tunneling rates are strongly suppressed as a function of s, as expected, and in 
general Unb > Uk- Moreover, in the lower panel, it is shown how inter particle interac- 
tions can be increased by considering stronger confinement along y—z directions. The 
regimes we will be mainly interested in, that is, the ones where magnetic properties 
may emerge, have to be characterized by sufficiently strong intra-species interactions 
and not too small hopping rates, so that robustness with respect to temperature is 
stronger. As such, intermediate values of s ~ 10 — 25 and possibly large values of sj_ 
will be our focus from here on. 



3 Perturbation theory at integer filling 

In this section, we analyse the magnetic phases deep inside the Mott regions at in- 
teger filling, njn, + tik € N, Uk = nm>, and discuss the qualitative phase diagrams 
realistically achievable with Rb-K mixtures. Our treatment flows along the lines of 
Ref. [15] , where the possible magnetic phases of ultra cold atomic mixtures have also 
been discussed. 

In the strong coupling limit (ti,t2 "C U^bK, U\, U%, where we have defined 1=K, 
2=Rb for compactness), processes in which single particle tunneling changes the total 
on-site populations require high energy and thus the low energy Hilbert subspace A 
contains only states with a fixed integer occupation number on every site. To derive an 
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Hamiltonian parameters 

tRb, tx, Uk, URb 

10- 7 U 




Fig. 1. Typical Hubbard parameters for 87 Rb- 41 K mixtures. Upper panel: Strength of the 
different Hamiltonian parameters in K units as a function of the lattice depth along the 
wire: tub (thin, red), £k (blue, dotted), Uk (dot-dashed, orange), and Um> (dashed, green). 
The tunneling rates decrease exponentially when increasing s. Lower panel: intra-species 
interaction scaling a function of the orthogonal lattice depth s± for different values of s — 
10, 15, 20, 25 from bottom to top. 



effective Hamiltonian acting in such subspace we use a generalization of the Schricffcr- 
Wolff transformation [35||36| . 

This kind of technique may be used whenever one has to deal with a Hamiltonian 
of the type: H — H +V, where the unperturbed Hamiltonian H may be diagonalized 
within a Hilbert subspace A [9j. In our case, V contains the hopping terms, whereas 
two-body interactions determine Hq. Here, A is the subspace with fixed occupation 
number for every site rij = v. By denoting with P, P 1 - the orthogonal projectors on 
the subspaces A, A 1 - respectively, one may write the Hamiltonian as the sum of a 
diagonal and an off-diagonal term: H = Ho + Hi with 

Ho = PH P + P^HoP 1 ^ (4) 
"Hi = PVP- 1 - +P ± VP + P ± VP ± . (5) 

Then one looks for a unitary transformation U = exp(ie<S'), with — S, that diag- 
onalizes the total Hamiltonian. The operator S may be found by imposing [Ho, S] + 
iH\ = 0, meaning that, up to the first non-trivial correction in e, the diagonal form 
of H is finally given by: 

H eS = H Q +i{S,H x }. (6) 



3.1 Spin-1/2 chains for v = 1 



Deep inside the first Mott lobe, when v = 1, the Hilbert subspace A is the tensor 
product of single-site Hilbert spaces which are two-dimensional and generated by the 
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two states: 

|ti) = atjvac) , \U) =a[ 2 \vac) . (7) 
On this space one can define the spin- 1/2 operators: 

■S? = (lti>Ui| + |4i)<ti|)/2 (8) 
S? = (|ti)<ii|-|W(ti|)/2< (9) 

s? = (\t i )(U-\U){U\)/2 (10) 

which allow to rewrite the Schrieffer- Wolff transformed Hamiltonian as: 

H = [J*S?S? +1 - Jj_ (SfSf +1 + SfSf +1 ) (11) 

i 

- ft (S? + Sf +1 ) - XI] 

where 

Jz = 2 tj + 4 < *4 J± = i^ ( 12 ) 

Un,bK U\ U2 UubK 

_ 2t| t 1 f 2 , + ^ /-] 

h -Th-u*' X -u[ + W 2 + w^- (13) 

Since in any canonical ensemble, /i is effectively fixed to (the effects of the trap- 
ping potential will be discussed below), we see that, up to a constant factor, the 



Hamiltonian represents a one-dimensional spin-1/2 XXZ model, as first found in 15 
Interestingly, this basic model displays already a rich phase diagram, characterized by 
phase transitions from a critical region for —1 < A = J z /J±_ < 1 to massive phases, 
either ferromagnetic (FP) for A < — 1 or anti-ferromagnetic (AFP) for A > 1. In the 
bosonic language, the former is a super-counter-flow (SCF) superfluid phase, which 
can be interpreted as a superfluid of particle-hole bound states, as easily seen by 
performing a particle-hole transformation in one of the species [14] . 

As shown in Fig. [2j the entire phase diagram can be scanned by just changing 
the interspecies interaction via Feshbach resonances. For very large values of UnbK, 
the system displays ferromagnetic behavior, corresponding in the atomic language 
to real-space phase separation; on the other hand, weak repulsion will indeed favor 
short-range antiferromagnetic interactions, due to the positive value of the diagonal 



interaction J z , as can be noticed from Eq. (13) 



3.2 Spin-1 chains for v = 2 and the Haldane phase 

In this section, we analyse the magnetic phases deep inside the Mott region with 
v = 2. In this case the low-energy Hilbert subspace is the tensor product of single-site 
Hilbert spaces which are generated by three states: 

(all) 2 (<4 2 ) 2 

|-i*) = ~^~l vac )' l +li ) = ^7^~l vac )' 

|0*> = 4i4 2 |vac) (14) 
which represent a multiplet for the spin-1 operators (15] : 

S+ = V5(|+li)(0 j | + |0 i )<-l i |) (15) 
Sr=y2(|-l i )<0 i | + |0 i )<+l i |) (16) 
Sf = |+l i )(+i i |-|_l i )<_i i | (17) 
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Fig. 2. v = 1,S = 1/2 case. Evolution of A — J z /J± with respect to the interspecies 
interaction U Rb n for s = 9, s± —65 (red, dashed line) and s = 15, s± — 65 (blue, solid 
line). The system can be driven from a ferromagnetic to a critical and, eventually, to an 
anti-ferromagnetic phase, by decreasing U Rb K- 



By means of the Schrieffer- Wolff transformation, one obtains an effective spin-1 
Hamiltonian which, neglecting constant and effective magnetic field terms, is given 
by: 



H cS = ^ v (Sf Sf +1 + SfSf +1 ) + \S*S* +1 + D(S, 

;0S?(fi 



z^2 

+ 5(^f +1 ) 2 + asf(s? +1 ) 2 + (sn 2 sui) (18) 



where — SfS? + S^Sf. In terms of the original parameters, the couplings are 
given by: 



V 



1 / tit 2 t\t 2 iifa tl~t 2 \ 2t\t 2 



2\U 1 ~2U RbK U 2 -2U RbK Ui U 2 ) U RbK 
3£? 3t§ t\ t 2 



X Ui U 2 + 2U RbK -U-l + 2U RbK - U 2 
D ^ Ui+U 2 -2U RbK t\ ( 3 



2 2 \U RbK -2U 1 Ui. U E 

t\ f 3 _8_ _ 5 

+ 2 \U RbK -2U 2 + U 2 U RbK 



3 + 1 7_ t 5 

U Rb K — 2U\ Ui — 2U Rb K Ui U Rb K 



2/3 1 75 

j £~ f _j_ _|_ 

\UnbK — 227 2 U2 — 2URbK U2 URbK 



ti 
+ i 



X =; 



\U Rb K — 
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^c/ 2 


— 2U Rb K 


3 
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U Rb K ~ 


t/i 
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U Rb K J 
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+ lV 


~ 3 l 
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-C/ 2 


U Rb K J 


, *1*2 
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t\t 2 t\t 2 t\t 2 t\t 2 2t\t 2 



2Ui 2U 2 2U x -AU RbK 2U 2 ~W RbK U RbK 

In Fig. [3j we show how the different parameters (measured in units of |7y|) depend 
on the interspecies interaction Uj^k for different values of s and s^. In a relatively 
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Fig. 3. v = 2, S = 1 case. The parameters of the effective Hamiltonian ( IS I as a function 
of the interspecies interaction Vbmk for s — 9, s± = 65 and assuming \rj\ = 1: A (black, 
dashed-dotted-dotted), D (orange, solid), 8 (green, dashed), ( (red, dotted) and \ (blue, 
dashed-dotted). Here, the hopping parameter is tx ~ 22nK. 



large parameter (Ur^r) region, <5, £, x are negligible with respect to A and D, so the 
Hamiltonian reduces to that of the so-called A — D model [3S|[42] : 

H X - D = ]T v(Sf Sf +1 + S?Sy +1 ) + \S?S* +1 + D{Slf. (19) 

i 

For example, if s = 15, s± — 65, this interval is given by 43nK < UnbK ^5 70nK, while 
if s = 9, sj_ = 65 it is given by 36nK < URbK ^ 60nK. 

The ground state phase diagram of the X — D model is even richer than the 
5=1/2 Heisenberg model one. Some qualitative considerations can be easily made 
when the first term in eq. |l9| ) can be neglected: if A is larger than D, the ground 
state reduces to a ferromagnetic state (if A < 0) or antiferromagnetic state (if A > 0); 
if D 3> |A| (large-D phase), spins in all sites prefer to have spin zero. For values of 



|D| which are comparable with |A|, instead, the first term in eq. (19) is important 
and new phases emerge, as sh own numerically in |37| , separated by critical lines with 
different universality classes [43]. More specifically^ if A is negative, the system is 
critical, belonging to the XX universality class, while for A positive a new massive 



phase arises, typical of isotropic spin-1 chain: the Haldane phase 44 . 

To characterize all massive phases, one may exploit the symmetries of the Hamil- 
tonian. Firstly, it has an explicit ^-symmetry to describe which one can consider the 
ferromagnetic order parameter and the Neel order parameter ((•) means expectation 
value in the ground state): 

<D$ M = lim (S?S?)-, (20) 
0£eci= I™ ((-l^'SfS?). (21) 

|j-j|->-oo 

As shown in [38] , the A — Z3-Hamiltonian has also a hidden Z 2 ^-symmetry which 
corresponds to a set of non-local order parameters, the so called string order param- 
eters [45] : 

O s Q tn, g = , I™ (-S?e^=^ s ?S?), (22) 

where a = x, y, z. The magnetic phase diagram can then be schematically represented 
as follows: i) in the ferromagnetic phase, only Of M is non zero, while all the other 
parameters are zero; ii) in the Neel phase, C^ ring = CSoci = for a = y, Op M = 
for all (3, while C| tring and 0^ 6cl are non zero; Hi) in the large-D phase, all the order 
parameters are zero. Finally, iv) in the Haldane phase, all three of 0" tring are non 
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Fig. 4. Phase diagram of Hamiltonian 



65 (upper) and s = 9, s± 



65 (lower panel). The dots represent the ground state of (19 1 by varying the interspecies 
interaction U^k (43nK< UnbK <70nK, lower 38nK< Ur^k <60nK) . When Urkk increases 
the ratio X/\r/\ becomes more negative in both cases. In both insets the parameter -q, fixing 
the energy scale, is plotted versus UmK- 



zero, while all of O^cci ano - °f ^fm are zer °j so that this phase can be seen as a 
dilute antiferromagnet, where antiferromagnetic correlations survive despite loosing 
the associated positional order present in the Neel phase. 

From the formulas given above, one can estimate the values of A and D and thus 
find that systems with s large (s > 15) will be mainly in the large-D phase, as it is 
shown in Fig. [4] (upper panel) for the case with s = 15, sj_ = 65 and UjibK in the 
interval ~ 43/70nK. When the value of s decreases, it is possible to go beyond the 
large-D phase and, in principle, reach the Haldane phase, as one can see from Fig. [4] 
(lower panel) for s = 9, s± = 65 and in a small interval of UubK- 

Nevertheless, we notice that the range of U^k in which the model reduces to 
the A — D one is not wide. Also, the Haldane phase region can be reached only for 
small values of s, where other terms in the effective Hamiltonian may compete with 
A, D and affect (or even destroy) topological order. Clearly, both these constraints 
represent a serious challenge for experimental realizations of the Haldane phase, as 
a clear observation would require both fine tuning of the parameters and very low 
temperatures. Finally, since the optical lattice depth in these cases would be of order 
s ~ 9, it would be preferable to check second order perturbation theory results against 
more accurate methods such as DMRG or Quantum Monte Carlo simulations. 
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Phase 


Ac 


As 


AFM 






SCF 







SF a + SF b 









Table 1. Strong coupling phases and relation with the charge and spin gaps. 



4 Numerical results 

In order to complement perturbation theory predictions, we will present here numer- 
ical results for the realistic Hubbard parameters discussed in Sec. [2] and equal filling 
for both species rim = uk — 1/2. 

Here, we are interested in i) determine a quantitative phase diagram spanning 
a broad range of interaction parameters, and ii) determine the optimal parameter 
regime in terms of resilience to thermal fluctuations, that is, identify the interaction 
strength around which the smallest between the charge and spin gaps, A c , A s is 
larger in absolute (that is, when expressed in nK units) value. This last observation is 
particularly interesting in view of possible experimental investigations, as it provides a 
quantitative guide to understand the critical temperatures needed to observe magnetic 
phases. Moreover, such information is not accessible within perturbation theory, as 
the largest gap region usually lies outside of its applicability regime, as shown below. 

All simulations have been carried employing the DMRG algorithm, a state-of-the- 
art method to tackle ground state properties of ID systems. In order to determine 
the quantitative phase diagram of the system and, at the same time, provide relevant 
energy scales for both charge and spin degrees of freedom, we have investigated the 
behavior of the charge: 

ML) = 1 \E L (L/2 + 1, L/2 + 1) + E L (L/2 - 1, L/2 - 1) - 2E L (L/2, L/2)] , (23) 

and spin 

A S {L) = X - [E L (L/2 - 1, L/2 + 1) + E L {L/2 + 1, L/2 - 1) - 2E L (L/2, L/2)} , (24) 

gap, where El(Ni, N2) is the ground state energy at size L in the particle number 
sector (TVi, AT 2 ). Far from the ferromagnetic regime UabK > URb, Ur, the extrapolated 
value of the two gaps A a = lim^oo A a (L) is sufficient to determine the magnetic 
phase diagram of the system according to Table [l] In particular, the Neel phase is 
fully gapped, while the super-counter-flow phase, correspondent to the critical phase 
of the XXZ model, has a finite A c but gapless spin excitations. 

In order to suppress boundary effects, we have employed periodic boundary condi- 
tions (PBCs); despite being computationally more requiring and limited to relatively 
small system sizes (here, up to L = 32), they assure much better scaling properties 
since translational symmetry is not explicitly broken, as reflected in the accurate 
finite-size scaling procedure (illustrated in Fig.|5|. High accuracy in ground state and 
excited state energies (with discarded weights usually of order 5 x or smaller) 

have been achieved by employing up to 1024 states per block and 5 finite-size sweeps, 
and further enhanced by applying an iterative procedure which used the final finite- 
size step at size L — 2 as the initial one at size L. Moreover, in each ground state 
sector, four states were usually targeted. In most simulation, a single-site basis of 8 
states (considering Rb particles as hard core) is sufficient to faithfully represent a 
realistic system; this feature has been systematically verified for small system sizes 
(L < 12), and by using a 15 state single-site basis close to critical points. 
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Fig. 5. Finite-size scaling analysis for both charge (upper panel) and spin (lower panel) 
gaps; here, s = 26, s± — 70, and inter particle interactions are varied. In both panels, lines 
are linear (Uri,k = 0.3) and algebraic (others) best fits. 



As a case study, we considered fixed values of the optical lattice depth, and in- 
vestigate the quantum phase diagram as a function of the inter-species interaction 
strength. This choice is motivated by the experimentally demonstrated tunability of 



Rb-K interactions 33 , which allows to span large interaction regimes by keeping 
the same lattice structure. In Table [2j we list the optical lattice setups investigated 
numerically, together with the corresponding hopping ratio. 

As already mentioned, PBCs assure very good scaling properties during the finite 
size procedure, where we have always considered polynomial scaling of the form: 

A a (l/L = x) = a Q + aix + a^a; 2 . (25) 

Different scaling ansatze with higher order contributions do not lead to significant 
quantitative changes. Typical results are shown in Fig. [5] illustrating how both gaps 
opens for relatively small coupling strengths. 

For all values listed in Table [2] we observe the expected phase diagram, which, 
as a function of UnbK, presents a fully gapless phase first, then an antiferromagnetic 
insulator, which melts into a SCF phase for relatively large values of the interaction 
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s 


20 


20 


26 


26 


Si 


50 


80 


70 


90 


t K (nK) 


5.4 


5.4 


2.8 


2.8 


tm/tK 


0.050 


0.050 


0.033 


0.033 



Table 2. Optical lattice configurations investigated in this section, together with correspon- 
dent hopping asymmetry and typical energy scales. 
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Fig. 6. Charge and spin gaps as a function of the inter particle interaction. While the charge 
gap increases significantly and then saturates as a function of Uri,k, the spin gap displays 
a non-monotonic behavior; it open almost simultaneously as A c (making the SCF region 
at small U very thin), displays a maximum around Um,K — 2, and then decreases till the 
BKT point separating the Neel phase from another SCF region. Numerical results show a 
good agreement with second order perturbation theory estimate (based on Eq. (12 1, ( 13 )) , 
indicated by a green arrow. 



strength. It is worth noticing that our procedure is not accurate enough to resolve 
a possible weak coupling SCF phase between fully gapless and AF phases; since the 
gap is expected to open according to Berezinskii-Kosterlitz-Thouless scaling [7], re- 
fined methods should in principle be employed in order to quantitatively tackle this 
issue. Nonetheless, such SCF phase will occupy a very small parameter region with 
extremely small charge gap, making its experimental observation extremely challeng- 
ing. 

Close to the SCF-AF critical point, we find excellent agreement between numerical 
and perturbation theory results, as illustrated in Fig.[6j we can then faithfully estimate 
the extent of the SCF region at strong coupling by using the results presented in Sec.[3j 
Depending on the values of s, sj_, the SCF is present in a region which is relatively 
wide, of order VbjokI^k — 2, making its strong coupling observation accessible in 
tunable systems. 

A crucial point in the search for magnetic phases in ultra cold atomic systems is the 
competition between thermal fluctuations and Neel order. At a qualitative level, one 
expects that larger spin gaps could provide stronger finite-temperature signatures of 
ground state physics, since thermal fluctuations would then be not sufficiently strong 
to activate spin excitations. In Fig. [7j we plot the spin gap in temperature units 
(K) for several choices of the lattice configuration as a function of the inter particle 
interactions. 
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Fig. 7. Charge and spin gaps (expressed in K units) as a function of the inter particle 
interaction. The maximum value of the spin gap is always reached far from the perturbative 
regime, making perturbation theory not fully reliable when estimating critical energy scales. 



Three main features emerge from the quantitative analysis. The first one is, very 
deep lattices favor larger values of A s , as expected, since they increase the ration 
between intra-species interactions and tunneling rates. Secondly, we observe a quasi- 
plateau structure of the spin gap, that is, its value is almost constant in a finite region 
of the AF phase; this implies that inter-species interactions do not need accurate finite 
tuning in order to find the optimal experimental configuration. Finally, it is worth 
noticing how the maximal value of A s is usually far from the perturbative limit, 
signaling how strong coupling perturbation theory is usually not a fully trustable 
tool in order to identify optimal experimental regimes where to search for AF. 

Finally, let us comment on the v — 2 case. As suggested by the results discussed 
in the previous section, the topological Haldane phase may emerge, but in very tiny 
parameter regimes; as such, it'd be interesting to investigate if such regimes, at least 
in the perturbative limit and eventually via more accurate numerical techniques, may 
enlarge in presence of larger mass imbalance by considering, e.g., 6 Li- 133 Cs mixtures. 



5 Discussion and conclusions 



We discussed magnetic phases in commensurate bosonic ultracold atomic gases, fo- 
cusing on the recently realized 87 Rb- 41 K mixtures in ID optical lattices. After an 
analysis on the relation between microscopic and Hubbard parameters, we presented 
a combined analytical and numerical study on the emergent insulating phases in such 
systems. By means of perturbation theory, we illustrated the general magnetic sce- 
narios accessible in such setups by comparing effective spin-chain Hamiltonians with 
the accessible parameter regimes, focusing on the relevant cases of S = 1/2 and 5=1 
spin representations. In the latter case, we showed how, while the so called large-D 
and ferromagnetic phases (phase separation) are well within reach, exploring topolog- 
ical states of matter such as the Haldane phase will represent a notable experimental 
challenge. In the former, we discussed how the entire strong coupling phase diagram, 
including ferromagnetic, antifcrromagnetic, and super-counter-flow phases, may be 
experimentally accessed by tuning the interspecies repulsion by means of, e.g., a Fes- 



hbach resonance 33 . As a quantitative benchmark for experimental realizations, we 
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have investigated the phase diagram of these systems by means of DMRG simula- 
tions. By evaluating the typical energy gaps in the systems, we estimated the optimal 
optical lattice setups in order to observe anti-ferromagnctism in cold atomic gases, 
finding that experimentally challenging but accessible temperatures 34 in the nK 
range are required to unambiguously reach such regimes. 
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